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Abstract 

In this paper the maximum likelihood equations for the parame¬ 
ters of the Weight Lindley distribution are studied considering differ¬ 
ent types of censoring, such as, type I, type II and random censoring 
mechanism. A numerical simulation study is perform to evaluate the 
maximum likelihood estimates. The proposed methodology is illus¬ 
trated in a real data set. 

Keywords: Weight Lindley distribution, Maximum Likelihood Esti¬ 
mation, Censored Data, Random Censoring. 


1 Introduction 

Advances in computational methods and numerical simulations have allowed 
to incorporate efficient models that are capable of describing real problems. 
Introduced by Ghitany et. al (2011) the Weight Lindley distribution with 
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two parameters, is very flexible model to be fitted by reliability data since 
this distribution has increasing and bathtub hazard shape. 

Some properties of this model were studied by Ghitany et. al (2011) as 
well as the parameter estimation based on the maximum likelihood method. 
Mazucheli et al. (2013) compare the efficiency of four estimation meth¬ 
ods: maximum likelihood, method of moments, ordinary least-squares, and 
weighted least-squares and conclude that the weighted least-squares method 
reproduces similar results to those obtained using the maximum likelihood. 
Using a Bayesian approach Ali (2013) consider different non-informative and 
informative prior for the parameters of the WL distribution. 

However, in studies involving a temporal response, is common the pres¬ 
ence of incomplete or partial data, the so called censored data (Lawless, 
2002). It is important to point out that even incomplete these data provide 
important information about the lifetime of the components and the omis¬ 
sion of those can result in biased conclusions. In literature there are different 
mechanisms of censorship (Balakrishnan & Aggarwala, 2000; Lawless, 2002; 
Balakrishnan & Kundu, 2013). Due to the large number of applications in 
medical survival analysis and industrial life testing, it will be considered cen¬ 
sored data with type II, type I and random censoring mechanism. Some 
referred papers regarding the reliability applications with those types of cen¬ 
soring can be seen in Ghitany & Al-Awadhi (2002), Goodman et. al. (2006), 
Joarder et. al. (2011), Iliopoulos & Balakrishnan (2011), Arbuckle et. al. 
(2014). 

The main objective of this paper is to estimate the parameters of the 
Weight Lindley distribution using the maximum likelihood estimation and 
considering data with different types of censoring, such as, type II, type I 
and random censoring mechanism. The originality of this study comes from 
the fact that, for the Weight Lindley distribution, there has been no previous 
work considering data with censoring mechanisms. 

The paper is organized as follows. In Section 2, we review some properties 
of the Weight Lindley distribution. In Section 3, we present the maximum 
likelihood method and its properties. In Section 4, we carry out inference 
for this model considering different censoring mechanism. In Section 5 a 
simulation study is presented. In Section 6 the methodology is illustrated in 
a real data set. Some final comments are presented in Section 7. 
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2 Weight Lindley Distribution 


Let T be a random variable representing a lifetime data, with Weighted Lind¬ 
ley distribution and denoted by WL(A, 0), the probability density function 
(p.d.f) is given by 


mx ^ )= (Mwm t(1+t)e ’ (1) 

for all x > 0 ,(ft > 0 and A > 0 and T(0) = J 0 °° e~ x x^ l dx is known as gamma 
function. The WL (jl| distribution can be expressed as a two-component 
mixture 

f(t\\,(f>) =pfi(t\X,4>) + (1 - p)f 2 (t\X,(ft), 

where p = A/(A + 0) and fj(t\ A, (ft) has p.d.f Gannna(0+j — 1, A) distribution, 
for j = 1, 2. 

The mean and variance of the WL distribution can be easily computed 
by 

_ </>(A + 0 + 1) 2 _ (</> + 1)(A + (ft) 2 — A 2 . . 

A(A + (ft) ,(T _ A 2 (A + (ft) 2 ' U 

The survival function of T ~ WL(6,c) with the probability of an obser¬ 
vation does not fail until a specified time t is 


S(t |A,0) 


(\ + (ft)r{(ft,\t) + (\t)*e- M 
(A + 0)T(0) 


(3) 


where T(x,y) = w y ^ 1 e~ x dw is called upper incomplete gamma. 

The hazard function quantify the instantaneous risk of failure at a given 
time t. The hazard function of T is given by 


f(t\X,<ft) A^ +1 t^ 1 (l + f)e xt 

S(t\X,(f>) (A + (ft)T((ft, Xt) + (Af)^ e~ xt 


(4) 


Figure [T] gives examples from the shapes of the hazard function for dif¬ 
ferent values of (ft and A. 
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Figure 1: Hazard function shapes for WL distribution and considering dif¬ 
ferent values of 0 and A 

3 Maximum Likelihood Estimation 


Using the classical approach, the maximum likelihood estimators was chosen 
due to its asymptotic properties. Maximum likelihood estimators are ob¬ 
tained from maximizing the likelihood function (see, Casella e Berger, 2002). 
The likelihood function of 6 — ($i,..., 6 k) given t, is 

n 

L(6,t) = l[mB). (5) 

i =1 

For a model with k parameters, if the likelihood function is differentiable 
at 9i, the likelihood equations are obtained by solving the equation system 

^-log(L(0,t)) = 0,i = 1,2, (6) 

The solutions of (4) provide the maximum likelihood estimators. In many 
cases, numerical methods such as Newton-Rapshon are required to find the 
solution of the nonlinear system. 

The maximum likelihood estimators of 6 are biased for small sample 
sizes. For large samples they are not biased and asymptotically efficient. 
Such estimators, under some regularity conditions, have an asymptotically 
normal joint distribution given by, 

(0) ~ iV fe [(0), / _1 (0)] para n ->■ oo, (7) 
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where 1(6) is the Fisher information matrix, k x k and Ijj(6), is the Fisher 
information of 0 in i and j given by, 

/ r\ \ 2 

I ij (6)=E — log(/ ( x \6))j , i,j = 1,2,...,k. (8) 

In the presence of censored observations, usually, its not possible to com¬ 
pute the Fisher information matrix, an alternative is consider the observed 
information matrix, where the terms is given by 

Hij(O) = (^ g9 d g9 log(. f(x\0))^j , i,j = 1,2,.(9) 

For large samples, approximated confidence intervals can be constructed 
for the individuals parameters dj, with confidence coefficient 100(1 — 7 )%, 
through marginal distributions given by 

(§i) ~ N[(9i), Hr7(6)} para n ->■ cxd. (10) 

4 Censoring and Parameter Estimation 

In this section, we provide the maximum likelihood estimator for the two 
parameters of the Weight Lindley distribution considering type II, type I 
and random censored data. Other types of censoring such as progressive 
type II censoring (Balakrishnan & Aggarwala, 2000) and Hybrid censoring 
mechanism (Balakrishnan & Kundu, 2013) can also be obtained to WL dis¬ 
tribution. 

4.1 Type II Censoring 

Usually in industrial experiments, the study of some electronic components 
are finished after a fixed number of failures r, in this case n — r components 
will be censored. This mechanism of censoring is call type II, see Casela & 
Berger (2001) for more details, and its likelihood function is given by 

1 r 

71' *1 —r 

£(A,#) = —f (H) 

' ' i =1 
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where is the order statistic. 

Let Ti, • • ■ , T n be a random sample of WL distribution, that is, T ~ 
hLL(A,0). The likelihood function is given by, 


L(X,(p\t) 


n \ \ r{<t>+1) ((A + 0)T(0, \t {r) ) + (A t {r) Y e 


(n — r)! 


(A + 


-X 


( 12 ) 


x 




e -A ho. 


i —1 


The logarithm of the likelihood function (12) is given by, 


01 1) = log(l + t (i) ) + (0 - 1) lo g(0b) _ A “ n log ( r (0)) 


2=1 


2=1 


2=1 


+ log(?r!) — log ((n — r)!) + r(0 + 1) log(A) — nlog(0 + A) 

+ (n-r) log ((A + ^)r(0, A« (r) ) + (A. 

(13) 

From <9/(A, c|£)/<9A = 0 and <9/(A, c|2)/<9c = 0, we get the likelihood equa¬ 
tions, 


n r(4> + 1) 


A + 4> 


A 


+ X] Hi) 

2=1 


(?r - r) (r(^>, A t {r) )e At M - (t (r) + 1) (At (r) )^ 


((A + </>)r(0, A t( r ))e t(r) ) + (A t( r )) 


(14) 


(n - r) ^r(0, At (r) ) + (A + 0)^(0, At (r) ) + (At^)^ log(At (r) )e At M 

((A + 0)r(^>,At( r ))) + (At( r )) <ji e At W 


+ 


n 


A + 4> 


+ ™ 0 ( 0 ) - ^ l0g ^«)’ 


= —rlog(A) 


(15) 


2=1 


where T(A:,a:) = J°° 1 log(u>)e "’die = log(a:)r(A;, a;) + xT(3, k, x), and 




a: ) is known as Meijer G-function. 


0 , 0,...,0 

s - 1 ,- 1,...,-1 

The solutions provide the maximum likelihood estimators of </> and A. 
Numerical methods such as Newton-Rapshon are required to find the solution 
of the non-linear system. 
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4.2 Type I Censoring 


In the presence of type I censored data, a fixed time t c is predetermined at 
the end of the experiment. Consider n patients in a treatment and suppose 
that d < n died until the time t c , then n — d are alive and will be censored. 
The likelihood function for this case is given by 

n 

L(0,t) = Ylme) Si s{t c \0) n - d , (16) 

i= 1 


where d = y~)"' Si is a random variable and Si — I (t t < t c ) is an indicator 
function. 

Let Ti, ■ ■ ■ ,T n be a random sample of WL distribution, that is, T ~ 
WL(A, 0). The likelihood function is given by, 


L(A,0|t) 


A^ +1 ) (A + 0)T(0, \t c ) + (Xt c y e 


<t> 0 -xt c 


n—d 


(A + 0)«r(0)' 


X 


X 




i e 


—Xt, 


Si 


i =1 


The logarithm of the likelihood function (17) is given by, 


n n n 

/(A, 0| t) = Si log(l + ti) + (0 - 1) ^ h lo g(^) - A ^ 

i— 1 2=1 2—1 

+ (71 - d ) log ((A + A t c ) + (Al c /e~ Atc ) 

+ d(<p + 1) !og(A) - 71 log(0 + A) - 77. log (r(0)). 


(17) 


(18) 


From dl( A, c\t)/dX = 0 and dl( A, c|2)/<9c = 0, we get the likelihood equa¬ 
tions, 


n d(0+1) 




(n - d) (V(0, A t c )e xtc - (t c + 1) (A t c )^ 


A + 0 A ^ ((A + 0)T(0, Xt c )e~ xtc ) + (At c 


(19) 


n - d) (T(0, A2 c ) + (A + 0)T(0, At c ) + (AtJ^log(At c )e 
((A + 0)T(0, At c )) + (Xt c )^ e~ Mc 


— 


= -dlog(A) 


7 








n 


n 


( 20 ) 


A + (ft 


+ nift((ft) - ^bgfe). 


2= 1 


The solutions provide the maximum likelihood estimators of (ft and A. 


4.3 Random Censoring 

In medical survival analysis and industrial life testing, random censoring 
schemes has been receive special attention. Suppose that the ith component 
could experiment censoring in time C t , then the data set is ( U,Si ), were 
ti = min(Tj,0 and Si — I (T < Cj). This type of censoring have as special 
case type I and II censoring mechanism. The likelihood function for this case 
is given by 

n 

L(e,t) = Y[f(t i \ 8 ) s ‘S(t i \ey- 5 ‘. ( 21 ) 

2=1 

Let Ti, ■ ■ • ,T n be a random sample of WL distribution, that is, T ~ 
WL(X,<ft). The likelihood function considering data with random censoring 
is given by, 

\ d (< j >+ 1 ) " / , X 1 - 5 , 

i(A ' m = (ITWt w n ( (A + A(i) + < Ai ‘) e ) 

* + t i )e- M )' i, . 

The logarithm of the likelihood function ( J22| ) is given by, 

n n n 

/(A, 01 1) = Si log(l + ti) + (<ft-1)^2 $i log (ti) - A ^ hU 
2—1 2=1 2=1 

n 

+ ^(1 - Si) log ((A + 0)T(0, A ti) + (A ti) V At ‘) 

2=1 

+ d(cft + 1) log(A) - n log(0 + A) - n log (T(0)). 

From dl( A, c\t)/d\ = 0 and dl( A, c|2)/<9c = 0, we get the likelihood equa¬ 
tions, 

n d((ft + 1) y, ^ ^ — Ai *) e — (U + 1) (Aii)^ 

A + ( ^ A 61 * * *=1 ((A + 0)T(0, \ti)e~ xti ) + (AL)^ 

(24) 

8 


x 


( 22 ) 


(23) 







n 


d l°g(A) 


E 


(1 - Si) (r((f), A ti) + (A + 0)T(0, A ti) + (Atj)^ log(Atj)e Ati ) 
((A + 0)T(0, A ti)) + (A ti)^ e~ xti 


n 

+ 7 - 7-7 + nii{<p) - V 5, log (ti) (25) 

A + ^ tr 

The solutions provide the maximum likelihood estimators of 0 and A. 


5 Simulation Study 


In this section we develop a simulation study via Monte Carlo methods. 
The main goal of these simulations is to study the efficiency of the proposed 
method. The following procedure was adopted: 

1. Set the sample size n and the parameter values 0; 

2. Generate values of the WL(0, A) with size n; 

3. Using the values obtained in step 2, calculate the MLE 0 e A; 

4. Repeat the steps 2 and 3 N times; 

5. Using § and 0 , compute the mean relative estimates (MRE) YliLi < Hf L i 

the mean square errors (MSE) YliLi , the bias n — anc ^ 

95% coverage probability. 

It is expected that for this approach the MRE’s are closer to one with 
smaller MSE. The 95% coverage probability was computed for the confidence 
intervals. For a large number of experiments, using a confidence level of 95%, 
the frequencies of intervals that covered the true values of 6 should be close 
95%. The type II censored data were drawn setting the completed data r 
and n — r were censored. To generate type I and random censored data, we 
utilize the same methods used by Goodman et. al. (2006) and Bayoud (2012), 
using these approaches it is expected that the proportions of censoring E\p\ 
are approximately 0.2 and 0.4. 

The results were computed using the software R (R Core Development 
Team). The seed used to generate the random values was 2014. The chosen 
values to perform this procedure were 6 = ((0.5, 2), (3, 2)), N = 100,000 and 
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n = (5,10,25,50,100). The values of 9 were selected to allow the increas¬ 
ing and bathtub shape in the hazard function. The maximum likelihood 
estimates were computed using the log-likelihood functions (13), (18) and 
(23) and the package maxLik available in R to maximize such functions. 
The coverage probabilities were also calculated using the numeric observed 
information matrix obtained from the maxLik package results. 

Tables 1-6 shows the MRE’s, MSE’s, Bias and the coverage probability 
C with a confidence level equals to 95% from the estimates obtained using 
MLE for N simulated samples, considering different values of n, 20% and 
40% of censored data. 


Table 1: MRE, MSE, Bias and C estimates for N samples of sizes n = 
(5,10, 25, 50,100), with 20% and 40% of type If censored data. 



MRE 

MSE 

Bias 

Cg5% 

MRE 

MSE 

Bias 

TWo 

n 

r 


4> = 

0.5 



A = 

2 


5 

4 

1.845 

0.810 

0.422 

0.969 

2.382 

30.240 

2.764 

0.961 

10 

8 

1.399 

0.199 

0.199 

0.971 

1.677 

8.414 

1.354 

0.965 

25 

20 

1.131 

0.033 

0.065 

0.963 

1.217 

1.243 

0.435 

0.962 

50 

40 

1.061 

0.012 

0.030 

0.957 

1.099 

0.393 

0.198 

0.957 

100 

80 

1.029 

0.005 

0.015 

0.953 

1.048 

0.157 

0.095 

0.954 

5 

3 

1.950 

0.990 

0.475 

0.968 

2.920 

55.550 

3.840 

0.950 

10 

6 

1.501 

0.281 

0.250 

0.971 

2.081 

19.370 

2.162 

0.959 

25 

15 

1.177 

0.049 

0.089 

0.968 

1.383 

3.220 

0.767 

0.962 

50 

30 

1.081 

0.016 

0.040 

0.960 

1.168 

0.867 

0.335 

0.960 

100 

60 

1.039 

0.006 

0.019 

0.953 

1.079 

0.309 

0.158 

0.956 

n 

r 


c\> = 

= 3 



A = 

2 


5 

4 

1.462 

10.280 

1.385 

0.950 

1.489 

4.874 

0.978 

0.955 

10 

8 

1.336 

5.116 

1.007 

0.959 

1.339 

2.243 

0.678 

0.961 

25 

20 

1.171 

1.802 

0.514 

0.959 

1.169 

0.753 

0.338 

0.961 

50 

40 

1.088 

0.740 

0.265 

0.958 

1.086 

0.303 

0.172 

0.958 

100 

80 

1.043 

0.304 

0.129 

0.955 

1.042 

0.125 

0.084 

0.955 

5 

3 

1.465 

12.250 

1.395 

0.944 

1.554 

7.022 

1.108 

0.947 

10 

6 

1.350 

5.994 

1.050 

0.954 

1.393 

3.159 

0.785 

0.956 

25 

15 

1.208 

2.342 

0.623 

0.960 

1.224 

1.158 

0.447 

0.961 

50 

30 

1.117 

1.048 

0.351 

0.958 

1.124 

0.499 

0.247 

0.959 

100 

60 

1.058 

0.428 

0.175 

0.956 

1.062 

0.201 

0.123 

0.957 
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Table 2: MRE, MSE, Bias, C and E[p\ estimates for N samples of sizes 
n = (5,10, 25, 50,100), with 20% and 40% of type I censored data. 



MRE 

MSE 

Bias 

CWo 

MRE 

MSE 

Bias 

cWo 

E\p\ 

n 


= 

0.5 



A = 

2 



5 

1.740 

0.844 

0.370 

0.950 

2.084 

25.800 

2.168 

0.926 

0.202 

10 

1.301 

0.159 

0.151 

0.961 

1.420 

4.533 

0.841 

0.946 

0.199 

25 

1.097 

0.028 

0.049 

0.958 

1.134 

0.799 

0.269 

0.949 

0.199 

50 

1.046 

0.011 

0.023 

0.955 

1.063 

0.313 

0.127 

0.949 

0.199 

100 

1.022 

0.005 

0.011 

0.952 

1.031 

0.139 

0.062 

0.951 

0.199 

5 

1.592 

0.536 

0.296 

0.947 

2.157 

27.270 

2.314 

0.909 

0.390 

10 

1.336 

0.205 

0.168 

0.957 

1.556 

7.256 

1.111 

0.929 

0.400 

25 

1.115 

0.039 

0.058 

0.959 

1.190 

1.478 

0.380 

0.941 

0.401 

50 

1.053 

0.014 

0.027 

0.956 

1.089 

0.563 

0.177 

0.944 

0.401 

100 

1.026 

0.006 

0.013 

0.952 

1.043 

0.247 

0.087 

0.948 

0.401 

n 


0 = 

= 2 



A = 

3 



5 

1.645 

25.380 

1.935 

0.903 

1.658 

11.720 

1.315 

0.902 

0.209 

10 

1.307 

5.907 

0.922 

0.941 

1.296 

2.455 

0.591 

0.941 

0.205 

25 

1.135 

1.593 

0.406 

0.954 

1.126 

0.623 

0.251 

0.954 

0.201 

50 

1.068 

0.663 

0.205 

0.954 

1.063 

0.261 

0.127 

0.954 

0.201 

100 

1.034 

0.287 

0.101 

0.953 

1.031 

0.115 

0.063 

0.952 

0.201 

5 

1.471 

17.960 

1.413 

0.906 

1.526 

9.547 

1.053 

0.903 

0.404 

10 

1.278 

5.844 

0.835 

0.937 

1.283 

2.581 

0.566 

0.935 

0.405 

25 

1.153 

2.076 

0.458 

0.951 

1.150 

0.908 

0.300 

0.949 

0.401 

50 

1.083 

0.914 

0.250 

0.954 

1.081 

0.403 

0.162 

0.953 

0.400 

100 

1.042 

0.390 

0.125 

0.953 

1.041 

0.176 

0.082 

0.952 

0.400 


It can be observed from the results that the Bias decreases as n increases 
and also as expected and the values of MRE’s tend to 1, allowing to get good 
inferences for the parameters of the Weight Lindley. That is, the estimators 
are asymptotically unbiased for the parameters. Moreover, the MSE of all 
estimators of the parameters tend to zero for large n, i.e. all estimators are 
consistent for the parameters. It is also important to point out that the 
coverage probabilities (C) of the two parameters approach the nominal one 
of 0.95 as there is an increase of the size of n. 
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Table 3: MRE, MSE, Bias, C and E[p\ estimates for N samples of sizes 
n = (5,10, 25, 50,100), with 20% and 40% of random censored data. 



MRE 

MSE 

Bias 

Cg5% 

MRE 

MSE 

Bias 

Cg5% 

m 

n 


= 

0.5 



A = 

2 



5 

1.887 

1.436 

0.443 

0.951 

2.282 

32.600 

2.564 

0.931 

0.200 

10 

1.330 

0.188 

0.165 

0.961 

1.490 

5.651 

0.980 

0.951 

0.201 

25 

1.104 

0.028 

0.052 

0.958 

1.152 

0.862 

0.304 

0.953 

0.201 

50 

1.048 

0.011 

0.024 

0.955 

1.070 

0.317 

0.141 

0.952 

0.201 

100 

1.024 

0.005 

0.012 

0.952 

1.033 

0.136 

0.067 

0.951 

0.201 

5 

1.868 

1.431 

0.434 

0.940 

2.492 

45.140 

2.984 

0.904 

0.409 

10 

1.381 

0.295 

0.190 

0.955 

1.636 

10.050 

1.271 

0.927 

0.418 

25 

1.117 

0.037 

0.058 

0.959 

1.198 

1.529 

0.397 

0.943 

0.419 

50 

1.053 

0.013 

0.027 

0.956 

1.091 

0.554 

0.182 

0.947 

0.419 

100 

1.025 

0.006 

0.013 

0.953 

1.042 

0.239 

0.085 

0.949 

0.419 

n 


0 = 

= 2 



A = 

3 



5 

1.960 

43.730 

2.880 

0.885 

1.938 

18.630 

1.875 

0.890 

0.193 

10 

1.441 

10.330 

1.322 

0.936 

1.419 

4.027 

0.837 

0.937 

0.202 

25 

1.162 

1.813 

0.486 

0.955 

1.153 

0.708 

0.305 

0.957 

0.203 

50 

1.077 

0.657 

0.231 

0.956 

1.073 

0.259 

0.146 

0.956 

0.203 

100 

1.037 

0.274 

0.111 

0.952 

1.035 

0.109 

0.070 

0.953 

0.203 

5 

2.040 

55.450 

3.119 

0.861 

2.086 

26.760 

2.171 

0.862 

0.382 

10 

1.527 

16.540 

1.580 

0.916 

1.520 

6.972 

1.040 

0.916 

0.404 

25 

1.189 

2.501 

0.568 

0.951 

1.184 

1.053 

0.368 

0.952 

0.407 

50 

1.091 

0.869 

0.273 

0.956 

1.089 

0.373 

0.178 

0.955 

0.407 

100 

1.043 

0.353 

0.130 

0.953 

1.042 

0.153 

0.084 

0.952 

0.407 


6 Application 


In this section, we illustrate our proposed methodology by considering two 
dataset. We will consider the Weight Lindley distribution to analyze such 
data. For sake of comparison we obtain the results with the Weibull and 
Gamma distributions and nonparametric KaplanMeier estimator (Kaplan & 
Meier, 1958). 

Firstly, to verify the behavior of the empirical hazard function it will 
be considered the TTT-plot (total time on test). Developed by Barlow and 
Campo (1975) the TTT-plot is achieve through plot of the values [r/n, G(r /n)\ 
where G{r/n ) = (£- =1 U + (n - r)t (r) ) /Ya=i U, r = 1 = 1 ,..., n 

and tM is the order statistics. If the curve is concave (convex), the hazard 
function is increasing (decreasing). When it starts convex and then concave 
(concave and then convex) the hazard function will have a bathtub (inverse 
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bathtub) shape. 

We also consider the AIC (Akaike Information Criteria) discrimination 
criterion method. Proposed by Akaike (1974) this method is based on the 
Kullback-Leibler Information. Let k be the number of parameters to be fitted 
and 6 the MLE’s of 0 , the AIC is computed through AIC = —2 log(L(0; t)) + 
2k. Given a set of candidate models for t, the preferred model is the one 
witch provide the minimum AIC value. 

6.1 Rats with vaginal cancer 

Presented by Pike (1966) the dataset is related to the lifetimes of 40 rats with 
vaginal cancer exposed to the carcinogen DMBA. In Table 4, we reproduce 
the data represented by survival times (in days) of 40 rats (+ indicates the 
presence of censorship). 


Table 4: Dataset related to the lifetimes of 40 rats with vaginal cancer ex¬ 
posed to the carcinogen DMBA. 


143 

164 

188 

188 

190 

192 

206 

209 

213 

216 

220 

227 

230 

234 

246 

265 

304+ 

216+ 

244 

142 

156 

173 

198 

205 

232 

232 

233 

233 

233 

233 

239 

240 

261 

280 

280 

296 

296 

323 

204+ 

344+ 


Based on Table 4, the data clearly has random censoring mechanism, 
consequently the equations (24) and (25) will be used to compute the MLE’s. 
Table 5 displays the MLE’s, standard-error and 95% confidence intervals for 
0 and A. Table 6 presents the results of AIC criteria, for different probability 
distributions. 


Table 5: MLE, Standard-error and 95% confidence intervals for 0 and A 


e 

EMV 

Erro P. 

IC 95 % (9) 

0 

21.7545 

1.3254 

(19.1566; 24.3523) 

A 

0.0978 

0.0066 

(0.0848; 0.1109) 
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Table 6: Results of AIC criteria for different probability distributions con¬ 
sidering the lifetimes of 40 rats with cancer. 


Criteria 

Weight Lindley 

Weibull 

Gamma 

AIC 

390.3428 

394.4234 

390.6482 


In the Figure [2j we have the TTT-plot, survival function adjusted by 
different distributions and KaplanMeier estimator and the hazard function 
adjusted by WL distribution. 





r/n 


Tempos 


Tempos 


Figure 2: TTT-plot, survival function adjusted by different distributions and 
KaplanMeier estimator and the hazard function adjusted by WL distribution 
considering the lifetimes of 40 rats with cancer. 

Based on the TTT-plot there is a indication that the hazard function 
has increasing shape. Comparing the empirical survival function with the 
adjusted distributions it can be observed a good fit for the Weight Lindley 
distribution. These result is confirmed from AIC since WL distribution has 
the minimum value. The hazard function adjusted by WL distribution con¬ 
firms the result obtained from TTT-plot. Therefore, through the proposed 
methodology the data related to rats with vaginal cancer can be described 
by the Weight Lindley distribution. 
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6.2 Lifetime of electrical devise 


Presented by Lawless (2002, p.112) the dataset is related to 60 electrical 
devices. The survival times is given in cycles to failure divided by 1000 
and was firstly presented without censoring. It will be considered that the 
experiment was ended after we observe r = 49 failure, therefore n — r — 11 
components will be censored. Table 7 reproduces the lifetimes from the first 
49 electrical devices. 

Table 7: Dataset related to the lifetimes of 60 (in cycles) electrical devices. 


0.014 

0.034 

0.059 

0.061 

0.069 

0.080 

0.123 

0.142 

0.165 

0.210 

0.381 

0.464 

0.479 

0.556 

0.574 

0.839 

0.917 

0.969 

0.991 

1.064 

1.088 

1.091 

1.174 

1.270 

1.275 

1.355 

1.397 

1.477 

1.578 

1.649 

1.702 

1.893 

1.932 

2.001 

2.161 

2.292 

2.326 

2.337 

2.628 

2.785 

2.811 

2.886 

2.993 

3.122 

3.248 

3.715 

3.79 

3.857 

3.912 







The experiment was ended after a predetermined number of failures 


therefore the data has type II censoring mechanism and the equations (24) 
and (25) will be used to compute the MLE’s. Table 8 displays the MLE’s, 
standard-error and 95% confidence intervals for cj) and A. Table 9 presents 
the results of AIC criteria, for different probability distributions, considering 
the electrical devices data. 


Table 8: MLE, Standard-error and 95% confidence intervals for cj) and A 
considering the electrical devices data 


e 

EMV 

Erro P. 

ICg 5 % (9) 

0 

0.6764 

0.1341 

(0.4137; 0.9392) 

A 

0.5260 

0.0954 

(0.3391; 0.7129) 


In the Figure [3j we have the TTT-plot, survival function adjusted by 
different distributions and KaplanMeier estimator and the hazard function 
adjusted by WL distribution. 

Similar to first dataset, comparing the empirical survival function with 
the adjusted distributions and through the AIC results it can be observed a 
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Table 9: Results of AIC criteria for different probability distributions con¬ 
sidering the electrical devices data. 


Criteria 

Weight Lindley 

Weibull 

Gamma 

AIC 

185.1739 

186.5965 

186.218 





r/n 


Tempos 


Tempos 


Figure 3: TTT-plot, survival function and the hazard function adjusted by 
Weight Lindley distribution considering the electrical devices data 


good fit for the Weight Lindley distribution. Based on the TTT-plot there is a 
indication that the hazard function has bathtub shape. The hazard function 
adjusted by Weight Lindley Distribution confirm those results. Therefore, 
through the proposed methodology the data considering the electrical devices 
can be described by the Weight Lindley distribution. 

7 Final Comments 


In this paper, we derived the maximum likelihood equations for the parame¬ 
ters of the Weight Lindley distribution considering different types of censor¬ 
ing, such as, type I, type II and random censoring mechanism. 

Based on simulation studies and on real applications, we demonstrated 
that using the proposed methodology it was possible to obtain good estimates 
of the parameters of Weight Lindley distribution. These results are of great 
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practical interest since this will enable us for the use of the Weight Lindley 
distribution in various application issues. 

There are a large number of possible extensions of the current work. The 
presence of covariates, as well as of long-term survivals, is very common in 
practice. Our approach should be investigate in both contexts. A possible 
approach is to consider the regression schemes adopted by Achcar & Louzada- 
Neto (1992) and Perdona & Louzada-Neto (2011), respectively. 
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